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^. Abstract 



We formulate a stochastic gauge fixing method to study the gauge dependence of the 
Abehan projection. We consider a gauge which interpolates between the maximal 
Abelian gauge and no gauge fixing. We have found that Abelian dominance for the 
heavy quark potential holds even in a gauge which is far from maximally Abelian one. 
The heavy quark potentials from monopole and photon contribution are calculated 
at several values of the gauge parameter, and the former part shows always the 
confinement behavior. 

1 Introduction 

Since 'tHooft and Mandelstam proposed the QCD vacuum state to behave like a mag- 
netic superconductor, a dual Meissner effect has been considered to play an essential 
role in the mechanism of color confinement ||l|, 0. A gauge is chosen to reduce the 
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gauge symmetry of a non-Abelian group to its maximal Abelian (MA) subgroup, and 
Abelian fields and magnetic monopole can be identified there. When one reduces 
SU{N) to U{1)^~^ by the partial gauge fixing, monopoles appear in ^(1)^~^ sector 
as a topological object. Confinement of QCD is conjectured to be due to condensation 
of the monopoles. By using the MA gauge which maximizes the functional 

R = ^Trf/^(x)a3?7^(a;)V3, (1) 

X,tJ. 

Suzuki and Yotsuyanagif^ first found that the value of Abelian string tension is close 
to that of the non-Abelian theory, where t/^(s) are SU{2) link variables on the lattice. 
Since then many numerical evidences have been collected to show the importance of 
monopoles in QCD vacuum: we refer to Ref.[0 for a review of these results. 

There are infinite ways of extracting U{1)^~^ from SU{N). This corresponds to 
the choice of gauge in Abelian projection. Abelian and monopole dominances can be 
clearly seen in MA gauge but not in the others; they seem to depend on the choice of 
gauge in the Abelian projection. However, the dual Meissner effect only in MA gauge 
is not enough for the proof of color confinement, since Abelian charge confinement 
and color confinement are different. 

Recently Ogilvie^] has developed a character expansion for Abelian and found 
that gauge fixing is unnecessary, i.e., Abelian projection yields string tensions of 
the underlying non-Abelian theory even without gauge fixing. Essentially the same 
mechanism was observed by Ambj0rn and Greensite for Z2 center projection of SU{2) 
link variables^. See also Ref.[^. Furthermore, by introducing a gauge fixing function 
Sgf = A X) TrU^(x)a"3U^(x) V3, Ogilvie has also shown that the Abelian dominance for 
the string tension occurs for small A. Hence he conjectures that Abelian dominance 
is gauge independent and that gauge fixing results in producing fat links for Wilson 
loop and is computationally advantageous for the measurements. Further Suzuki 
et.al. have shown that if the gauge independence of Abelian dominance is realized, 
the gauge independence of monopole dominance is also proved |§]. Hence to prove 
the gauge independence of Abelian and monopole dominances are very important 
especially in the intermediate region between no gauge fixing and exact MA gauge 
fixing. 

In this letter, we analyze the gauge dependence of the Abelian projection. We now 
employ stochastic quantization with gauge fixing as the gauge fixing scheme which 
has been proposed by ZwanzigerP] 

where x is Euclidean space-time and r is fictions time. 77 stands for Gaussian white 
noise 

K{x,t)) = 0, 



{r^l{x,T)ril{x'y)) = 25'^%j\x-x')5iT-T'). 
Here A is defined as 

A{x) = A\x)ai + A'^{x)a2, 

A\x) = -^gT.[d,K(^) - 9Alix)Alix)], (3) 

A\x) = -A9j:[d,K(^) + 9Al{x)Al{x)]. (4) 

Note that a = corresponds to the MA gauge fixing and a = cxd is the stochastic 
quantization without gauge fixing. 

2 Formulation 

SU{2) elements can be decomposed into diagonal and off-diagonal parts after Abelian 
projection, 

Uf,{x) = c^{x)Uf,{x), (5) 

where c^(x) is the off-diagonal part and u^{x) is the diagonal one, 

( \ — /^6xp(26'^(x)) \ 

""^^"^'-y exp{-ze,{x)) ) ■ 

The diagonal part can be regarded as link variable of the remaining U(l). One can 
construct monopole currents from field strength of U(l) links [p!0|: 



0^,uix) = 9^{x) + 9^{x + ij,)-e^{x + 0)-9^{x) 
= 6f,^ + 27m^^ (-TT < 6f,^ < tt), 

k^,{x) = -—e^^p^djp^{x). (6) 

Wilson loops from Abelian, monopole and photon contributions can be calculated 



as [y] 



^Abelian ^ exp(-| J] M^,(x)^^,(x)), (7) 

^monopole ^ ^^^^^^^ ^ kp{s)D{x - x')^e^(,p^d^Mp,{x')) , (8) 

x,x',a,l3,p,o- 

^photon ^ g^p^_^ ^ d;e,,ix)Dix-x')Mx')), (9) 

x,x',fj.,u 

Ju{x) = d~M^^{x), 
3 



where 9 is a lattice forward derivative, d~ is a backward derivative and D{x — x') is 
the lattice Coulomb propagator. J^ is the external source of electric charge and M^j, 
has values ±1 on the surface inside of Wilson loop. In order to achieve better signals, 
we perform smearing for spatial link variables of non-Abelian, Abelian, monopole and 
photon|jl2[ .We set 7 = 2.0 for non-Abelian configurations and 7 = 1.0 for the others, 
where 7 is the parameter which determines the mixing between a link variable itself 
and staples surrounding the link 

f/^(x)— 'i = ^{7 U,{x) + 5:(staples)}. 

Here A^ is a normalization factor. 

Stochastic quantization is based on Langevin equation which describes stochastic 



processes in terms of fictions time[|T^. A compact lattice version of this equation with 

(10) 



gauge fixing was proposed in Ref.|ll^ 

U^{x, T + 6t) 



J fl 



UJ(X,T) 



uj{x, Tyexp{if^aa)Uf,{x, t)uj{x + fi,T), 
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^L{x,r)aaST). 



In MA gauge. 



Aiat(a:,T) 



i[a3,X{x,T)] 

2(X2(x, r)cri - Xi{x, r)cr2) 



where 



X(x,t) 



^(f/^(x, r)a3U^{x, tY -U^{x- fi, r) ^3^/^(2; - /i, r)) 



(11) 

(12) 
(13) 



As an improved action to reduce finite lattice spacing effects, we adopt the Iwasaki 
action iH 



P 






S=2-^Tr(Co£| ^ t+Ci^t^ 



'0 



where Co + 8C1 = 1 and Ci = —0.331. The Runge-Kutta algorithm is employed for 
solving the discrete Langevin equation [|T6|. As will be shown, the systematic error 
which comes from finite St is much reduced. 



3 Numerical Results 

Numerical simulations were performed on 8^ x 12 and 16^ x 24 lattices with jS = 
0.995, a = 0.1,0.25,0.5,1.0, and 5t = 0.001,0.005,0.01. Measurements were done 
every 100-1000 Langevin time steps after 5000-50000 thermalization Langevin time 
steps. The numbers of Langevin time steps for the thermalization were determined 
by monitoring the functional R and Wilson loops. 

In Fig.|l], we plot (A^)^ + (A^)^ as a function of the gauge parameter a. a = 
corresponds to A = 0, i.e., MA gauge. When a increases, the deviation from the 
gauge fixed plane becomes larger. 
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Figure 1: Gauge parameter a versus (A^)^ + (A^)^. Lattice size is 8'^ x 12, 6t = 0.005 
and P = 0.995. 

We calculated the heavy quark potentials from non-Abelian, Abelian, monopole 
and photon contributions by 

V(R) = - lim log , <^(-^'^)) (14) 

We fit them with the following function, 

V{R) = Vo + aR+^. (15) 

ri 

In order to check that our Langevin update algorithm with the stochastic gauge 
fixing term works correctly, we plot in Fig.|] the heavy quark potential V{R), which 
is consistent with the result by the heatbath update. Runge-Kutta method works 
well, i.e., we see no difference among data with 5t = 0.001,0.005 and 0.01. 

In Fig.|^ we show the Abelian heavy quark potentials for different a's together 
with that of non-Abelian potential. The heavy quark potentials from monopole and 
photon contributions are plotted in Fig.|[ They can be well fitted by a linear and 
Coulomb terms, respectively. We see that the linear parts of potentials are essentially 
same from a = 0.1 to 1.0, and all of them show the confinement linear potential 
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Figure 2: Non-Abelian heavy quark potentials derived from heatbath and Langevin 
updates for three values of r. Lattice size is 8^ x 12 and (3 = 0.995. 
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Figure 3: Heavy quark potentials from non-Abelian and Abelian contributions. 
Lattice size is 16^ x 24, 6t = 0.005 and /3 = 0.995. 



behavior. Therefore even when we deviate from the MA gauge fixing condition, we 
can identify the monopole contribution of the heavy quark potential showing the 
confinement behavior. As a increases, statistical error becomes larger. This result 
suggests that the gauge fixing is favorable for decreasing numerical errors as pointed 
by Ogilvie||]. 

In Fig.|^ we plot the values of the string tensions from Abelian, monopole and 
photon contributions as a function of the gauge parameter a. They are obtained 
by fitting the data in the range 2.0 < R < 7.0. We have taken into account only 
statistical errors. The upper two lines stand for the range of the non-Abelian string 
tension. The Abelian and the monopole dominances are observed for all values of a. 
The string tensions from the Abelian parts are about 80% of the non-Abelian one. 
We expect that the difference of the percentage between our result and that of Bali 



et.al.[0 becomes smaller when we go to larger lattice size. On the other hands, the 
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Figure 4: Heavy quark potentials from monopole (left) and photon (right) contribu- 
tions. Lattice size is 16^ x 24, 6t = 0.005 and P = 0.995. 
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Figure 5: String tensions from non-Abelian, Abelian, monopole and photon contri- 
butions. Lattice size is 16^ x 24, 6t = 0.005 and P = 0.995. 

string tension from the photon part is consistent with zero. 

4 Concluding Remarks 

We have developed a stochastic gauge fixing method which interpolates between the 
MA gauge and no gauge fixing. In Refs.|]^,|0 and [0, effects of gauge fixing are 
studied for lattice algorithms where gauge fixing is done after field configurations are 
generated. In the stochastic gauge fixing procedure studied here, the attractive force 
to the gauge fixed plane along a gauge orbit is applied together with the Langevin 
update force. The method is tested together with the Iwasaki improved action and 
the Runge-Kutta algorithm. We have found it works well. 

We have studied the gauge dependence of Abelian projected heavy quark poten- 
tial. It is observed that the confinement force is essentially independent of the gauge 
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parameter. In the calculation of Abelian heavy quark potential, we have seen that as 
gauge parameter a increases, the statistical error becomes larger. This result suggests 
that the gauge fixing is favorable for increasing the statistics as pointed by Ogilvie0. 
It is expected that as a increase, Abelian string tension would approach the non- 
Abelian one p, ^. Therefore it is important to see behavior of the string tension 
as a becomes much larger than one. But data are more noisy for large a and we 
are planning to employ a noise reduction technique such as integral method |0 for 
obtaining statistically significant data. 

It is desirable to study the gauge dependence (or independence) of other quantities, 
such as the monopole condensation, which may reveal the role of gauge fixing in the 
dual superconductor scenario. Hioki et.al. have reported the correlation between 
monopole density and Gribov copy in MA gauge fixing [|1^]. Gribov copy effect for 



Abelian string tension has been studied by Bali et.al. |12|. For Landau gauge, Gribov 



copy may be avoided by Langevin update algorithm together with the stochastic 



gauge fixing|14]. It is, therefore, interesting to investigate effects of Gribov copy in 



Abelian projection with the method. 
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